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Many important transport phenomena are described by simple mathematical models rooted in 
the diffusion equation. Geometrical constraints present in such phenomena often have influence of 
a universal sort and manifest themselves in scaling relations and stable distribution functions. In 
this paper, I present a treatment of a random walk confined to a half-space using a number of 
different approaches: diffusion equations, lattice walks and path integrals. Potential generalizations 
are discussed critically. 
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I. INTRODUCTION 

Diffusion is perhaps one of the most well-studied fields 
of modern physics. Diffusion processes are ubiquitous 
and extremely important in physics, biology, chemistry, 
material science and other disciplines of basic and applied 
research. Many famous names, such as Fourier, Kelvin 
and Einstein, are associated with the theory of diffusion. 
About two hundred years after the diffusion equation was 
first written down, the motivation for studying diffusion 
by both theorists and experimentalists has not faded. 
Modern research deals with questions of much greater 
complexity than the originally formulated diffusion equa- 
tion. Fields such as pattern formation, non-equilibrium 
growth, reaction-diffusion processes and granular media 
keep challenging physicists and applied mathematicians 
and yield many interesting and often unexpected results. 

Compared to these actively developing fields, "old- 
fashioned" linear diffusion may seem boring and trivial. 
Often it appears though, that newly discovered phenom- 
ena and incessantly developing techniques provide formal 
and conceptual contexts in which old and well-known 
questions appear in a new light and acquire additional 
scientific and pedagogical value. For instance, rapid 
progress in the physics of polymers and critical phenom- 
ena in the late 1960's and 1970's, stimulated by introduc- 
tion of field-theoretic ideas and methods, has shifted the 
accents in diffusion studies towards scaling and universal- 
ity. Many systems that appeared completely unrelated 
at first glance turned out to exhibit identical behavior in 
certain asymptotic ( "critical" ) regions of the parameter 
space. 

In this paper, we study the rather old problem of a 
random walk 13 confined to a half-space, using a number 
of different approaches. First, we review the traditional 
methods - the diffusion equation and lattice walks. Next, 
we discuss the applications of these methods to Gaus- 
sian polymer chains, with emphasis on scaling. Finally, 
we show how the same problem can be solved using ba- 
sic field-theoretic methods, i.e. functional integration in 
Fourier space. We shall see how important scaling rela- 



tions and universal distribution functions can be obtained 
in each case through appropriate limiting procedures. We 
also discuss possible ways to extend and generalize the 
problem for more realistic systems. 



II. DIFFUSION AND RANDOM WALKS 

A. Boundary conditions 

Consider a -ZV-step random walk starting at ro = in 
the three-dimensional (3D) space. Let G(r, TV) be the 
probability density for the walk to end at r. It is well 
known that for large N and in the absence of obstacles 
and boundaries, G(r, N) is a solution of the diffusion 
equation 



(d N -V 2 )G(r,N) = 6(r)6(N), 



namely, 



G(t,N) 



e -r 2 /(2N) 

(2nN) 3 / 2 ' 



(1) 



(2) 



It is reasonable therefore to assume that solutions of the 
same kind can as well be found for any bounded region 
R. 

Naturally, one has to specify the boundary conditions. 
This is not as trivial as it appears and, in fact, depends on 
the physical context of the problem. If, for example, the 
random walker is allowed to touch the boundary and then 
step back with probability one, the reflecting boundary 
conditions are appropriate. Formally, it means that the 
flux across the boundary vanishes 



VG(r,iV)-n| rGaR = 0. 



(3) 



Here SR denotes the boundary of the region R and n 
is a unit vector locally normal to <9R. Another possible 
choice of boundary conditions corresponds to the case 
when the walker sticks to the boundary upon reaching it 
- the absorbing boundary conditions 



G(r,A0| re9R = 0. 



(4) 



For reasons that become clear below, we shall mainly be 
interested with the latter case. 



B. Method of images 



heat-conduction problems with a variety of sources and 
boundary conditions 3 . 

Kelvin was also the first to apply the method of im- 
ages to account for boundary conditions for electricity 
conduction in a semi-infinite telegraph line 2 . 



Consider a random walk starting at ro = za away from 
the plane z = and confined to the z > half-space. For 
the absorbing boundary, we expect the probability distri- 
bution for the end point to satisfy the following boundary 
value problem 



G(z 



V 2 )G(r,N;a) = S(r 
:0) = 0, G(|r|- 



3C 



za)S(N) 




(5) 



Any introductory textbook on PDEs contains a straight- 
forward solution of this problem, which consists of in- 
troducing a sink, or negative image, at (0, 0, —a) and ex- 
tending the problem to the entire space«. The boundary 
condition at z = is then automatically satisfied and 
the solution is 



G(r,N;a) 



1 



(2ttA0 3/2 



,-(r-i.a) 2 /(2N) 



-(r+ia) 2 /{2N) 



(6) 

For a <C \/N, we can expand the expression in parenthe- 
ses to obtain 



,'2m 
G(r,N;a)± ( — 



e -r 2 /(2JV) 

(2ttA0 3 / 2 ' 



(7) 



We see that the probability distribution can be factor- 
ized into z-dependent and z-independent parts. The 
latter, which includes degrees of freedom parallel to the 
boundary, is not affected by the presence of the boundary. 
Thus, in what follows we will be predominantly occupied 
with the z-dependent part of G(r, N; a). 



C. A short historical digression 

The theory of diffusion was first developed in the be- 
ginning of the 19th century by Joseph Fourier; his work 
was summarized in the famous Theorie analytique de la 
chalew^, first published in 1822. It contains an exten- 
sive treatment of homogenous heat diffusion problems for 
a variety of geometries, mostly by the variable separation 
method. 

The first generalized approach to solving non- 
homogenous diffusion probems was formulated by Sir 
William Thomson^, more widely known as Lord Kelvin, 
in 1850. He realized that particular solutions can be ob- 
tained by superposition of solutions for "instantaneous 
simple point sources" (which are now called by physi- 
cists "Dirac's delta-functions"). In short, what he did 
was to invent the Green's function method for the diffu- 
sion equation; it was later used by E. W. Hobson to treat 



D. "Phantom" polymers and scaling 

As mentioned in the Introduction, random walks are 
often used to model long polymer molecules. Proper de- 
scription includes a nontrivial requirement: the random 
walk must be self-avoiding, i.e. no point in space should 
be visited more than once. This requirement introduces 
long-range correlations and makes the problem tractable 
only approximately. However, if the self-avoidance is 
removed, the problem becomes much simpler. Such a 
polymer is called phantom; alternative names are ideal or 
Gaussian chain. The distribution function for the end- 
to-end radius vector r of a phantom polymer of length 
N obeys the diffusion equation Q. 

Rapid developments in polymer physics and the theory 
of critical phenomena have revealed a number of univer- 
sal properties that arise in all polymer chains beyond 
a certain level of coarse-graining 4,5 . These properties 
are characterized by scaling relations. Perhaps the most 
widely known scaling law relates the mean end-to-end 
distance R of a polymer chain to its length: R oc N". 
The number 1/v thus plays the role of the polymer's 
fractal dimension. It is universal in that it depends only 
on the dimensionality of the embedding space, e.g. in 
3D, v w 0.59 for a self-avoiding polymer and v = 1/2 
for a phantom one. Another important scaling relation 
describes the number of different configurations *Tt of a 
polymer 



m = const x C N 7 



(8) 



where C is the "effective coordination number" that de- 
pends on the microsopic details (cf . Eq. fH|) , and 7 is 
a universal exponent. The factor ( N can be thought 
of as counting the configurations of an unconstrained 
iV-step random walk with £ options available at each 
step, whereas iV 7-1 accounts for constraints such as self- 
avoidance, obstacles present etc. 

If the distribution function G(r, N) is known, the expo- 
nent 7 can be obtained in a very straightforward way, by 
simply integrating G(r, N) over the whole space. Thus, 
a phantom polymer has 7=1. It turns out that incorpo- 
rating self-avoidance constraint leads to 7 = 1.161*1. This 
can be interpreted as the enhancement of available space 
for a self-avoiding polymer which appears "swollen" com- 
pared to a phantom one. 

Consider a phantom polymer anchored to the xy plane 
and confined to the z > half-space. Assuming that the 
polymer is strongly repelled from the z = plane, what is 
the corresponding boundary value problem? Intuitively, 



it would seem that we must adopt reflecting (zero flux) 
boundary conditions at z — 0. However, a more careful 
analysis shows that this is incorrect. Note that G(r,N) 
is proportional to the number of ./V-step paths of length 
iV ending at r. Since the plane is repelling, while count- 
ing the paths contributing to G(r, N), we should discard 
those touching the z = plane at least once. Thus, 
the probability and not the probability flux should van- 
ish at the boundary, which corresponds to the absorb- 
ing boundary conditions. This is quite counter-intuitive: 
while a polymer, being an entire path, is repelled from 
the boundary, a fictitious random walker that we employ 
to model this path is absorbed there! 

Chandrasekha»& in his classic paper of 1943, considered 
both kinds of boundary conditions (we will touch on his 
derivation below). However, the abovementioned sub- 
tlety was overlooked by a number of subsequent authors 
which lead to incorrect calculation of average quantities 
in polymer adsorption problems (for a very lucid discus- 
sion of this topic see the paper by DiMarzioi) . 

Thus, the anchored Gaussian chain should be described 
by Eq. ©. Apparently, we should take a = + , which 
leads to 



G(r,iV;0+) = 0. 



(9) 



Of course, this result does not make sense. To obtain the 
correct result, we should have normalized the distribution 
function for a finite value of a and only then take the limit 
a — > + , which would yield 



/ z \ e~ r /( 2JV ) 
G(r,N;0+) = (-) — — . 



(10) 



However, normalization eliminated any information con- 
nected with 7; indeed, any normalized distribution func- 
tion by definition has J d 3 r G(r,N) = 1. Again, one 
can see that the problem here is that of the order of tak- 
ing limits. For any physical polymer a is finite, at least 
of the order of the smallest coarse-graining scale - the 
persistence length. For small but finite a, 



m<x 



dx dy 



dz 



N 5 / 2 



e-^ocaJV- 1 / 2 . (11) 



From here, we can read off the value of 7 which for the 
anchored random walk is denoted 71—: 



71 



(12) 



Thus, we see that the presence of the plane reduces the 
number of accessible configurations compared to the un- 
constrained case which manifests itself in the scaling ex- 
ponent 71 . Note that reflecting boundary conditions pro- 
duce 7 = 1 so that the number 91 is unchanged relative 
to the unconstrained case. This is yet another argument 
for incorrectness of reflecting boundary conditions in this 
problem. 



E. Counting walks on a lattice 

Chandrasekhar^ suggested a direct way of counting 
the paths on a lattice when a reflecting or an absorb- 
ing boundary is present. We will briefly describe the 
derivation for an absorbing boundary. The reader is en- 
couraged to read the original paper which despite being 
written more than half a century ago, remains one of 
the best introductions into random walks and stochastic 
processes in general. 

Consider a one-dimensional random walker on a lat- 
tice (discrete z-axis) with absorbing boundary at z = 0. 
Suppose, the walk starts some distance n from the origin; 
our task is to calculate the number of paths leading from 
n to some other point m, without touching the boundary. 
It turns out that it is easier to calculate the number of 
paths that do touch the boundary and then to subtract 
it from the total number of paths leading from n to m. 
To do so, we make use of a very elegant theorem - the 
reflection principle. 



-n 







m 



N 




FIG. 1: The reflection principle. 



Let us extend our lattice to include the negative part 
of the z-axis as well. Then, the reflection principle states 
that the number of iV-step paths originating at n, ending 
at m and touching or crossing the boundary z = is equal 
to the number of iV-step paths that originate at — n and 
end at m. Figure ^ illustrates the reflection principle by 
presenting a way to build a one-to-one mapping between 
the two sets of paths. Thus, the number of paths not 
touching the boundary is 



m = 



N 



l[N + m- 



N 



\[N + m + n] 



(13) 



For the starting point near the boundary, and m <C N, 
we can expand the binomial coefficients using Stirling's 



formula to obtain 

m~ 2 N 



2 \ 1/2 to e - m2 /( 2jy ) 



itN 



N 



(14) 



Dividing by the total number of paths of length N (which 
is 2^), we obtain the probability density for a path to 
start near the boundary and to end at some point m 
without returning to the boundary 



G(m, TV) 



2m e -m 2 /(2JV) 

~W (2TTN) 1 / 2 ' 



(15) 



III. RESULTS FROM PATH INTEGRALS 

In this section, we employ the methods of functional 
(path) integration to build a field-theoretical model that 
describes a Gaussian chain anchored to an impenetrable 
plane. 



To calculate this sum, one has to define a measure of 
integration. One way would be to discretize the chain 
and view Z(N) as a limit of a multidimensional inte- 
gral. In this case, the problem is almost identical to 
calculating a free particle propagator a Id Feynman & 
Hibbsiii. Another way around is to count the Fourier- 
components of c(t); this way is somewhat easier, since 
for harmonic Hamiltonians, degrees of freedom decouple 
in Fourier space. Thus, setting 



D[c(r)} 



n 



d 3 c(q) 
(2tt) 3 



(19) 



and Fourier-transforming Equation l|18fl . we obtain 



G(k,N) = 



TT / rf3£ (g) c »k.5(q)k-^-ll-lq 2 |cfq)| 2 

H / (2tt)3 

11 J (2tt)3 



(20) 



A. Example: unconstrained Gaussian chain 

In the field-theoretic approach, a flexible chain is de- 
scribed by a function ("path") c(r), where r measures 
the position along the chain. The energy of a self- 
avoiding chain in an external potential is given by^S: 



v 

+ 2 



JV 



JV 



c 2 (t) dr+ / U[c[r))d- 



N /-JV 



JO 



S[c(t) - c(r')]drdr'. 



(16) 



The first two terms can be viewed as a harmonic poten- 
tial between neighboring segments of the chain and the 
external potential, respectively, whereas the last one ac- 
counts for excluded volume effects: each time the chain 
self-intersects, a penalty of v is paid. In what follows, we 
omit the self-avoidance constraint. 

The partition function for such a chain is a sum over 
all possible paths c(r), given by a path integral 



Z(N)=JD[c(r)} 



,-H\, 



(17) 



If the external potential is set to zero, this expression 
simply counts the number of configurations of a Gaus- 
sian chain of length N. The configurations are weighted 
with a weight e~ H ^. If we want to count only paths 
starting at the origin and leading to some point r then, 
after normalization by Z(N), we obtain the probability 
density 



G(r,iV) 



1 



Z(N) 



D[c(t)} 6[c(N) 



c{0)-r}e- H[c] . 
(18) 



Both the numerator and the denominator contain prod- 
ucts of Gaussian integrals that can be calculated by 
"completing the square." The result is 



G(k, N) = exp k 



dq 1 — cos qN 



2tt q 2 
which is the Fourier-transform of 

e -r 2 /(2N) 



G(r,N) 



(2nN) 3 / 2 ' 



_ e -k 2 W/2 j 

(21) 
(22) 



B. Anchored polymer - the partition function 

Now that we are somewhat familiar with the method- 
ology of path integrals, we return to the original problem 
- the anchored Gaussian chain in a half-space z > 0. 
Writing it in terms of path integrals, we immediately see 
that our task will not be as simple as before. The reason 
for this is that now possible values of c 2 (t) should be pos- 
itive. If we stay in real space and calculate the limit of 
a multidimensional integral, we see that "completing the 
square" does not work because of this constraint, since 
the resulting integrals cannot be calculated analytically. 
If we decide to move to Fourier space, it is not even clear 
how to define the measure of integration. 

The suggested way out of this complication is as fol- 
lows. We allow the polymer to cross the boundary and in- 
troduce a strong repulsive interaction between the plane 
at z — and the chain. Each time the polymer crosses 
or touches the plane, it is "penalized" by a large amount 
of energy. The modified Hamiltonian is then 



H = H a + H 1 



(23) 



where 



and 



H °=2 



N 



N 



c 2 dr. 



H 1= g S[c z (t) - c z (0)} dr. 



(24) 



(25) 



Thus, we expect that when the coupling constant g > 
becomes infinitely large, the polymer will be entirely on 
one (either positive or negative) side of the plane. The 
partition function is then 

Z(g,N)= [ D[c(t)] g-Ho[c]-fl- 1 [ ff ,c] - (26) 



To evaluate Z(g,N), we expand the integrand in 
Eq. (|26JI in powers of g. Such an expansion could be 
problematic when g is large, and this is the limit we are 
primarily interested in. However, if we are able to cal- 
culate the general term of the expansion and to perform 
the summation to infinity, this approach is valid. 

The n-th (n = 1, 2, ...) term of the expansion reads 



(-.9)' 



/=i 



D[c(t)] e- H °M Y[ / 6{c z (n) - c z (0)] dn. 



Ordering the set {t;} and Fourier-transforming the <5- 
functions, this can be rewritten as 



2^] ID[c(r)]e- 



- Ho[c] / dr n 
Jo Jo 



i=i 



dr„_i... / dn f[ / dh e -**l°'M-°'<f»]. 



(28) 



Henceforth, we focus on the c z (t) and denote it c(t) for the sake of simplicity. Integrating it out by the method 
discussed above (i.e. "completing the square"), we are left with 



n r N 



— ) / dr n J dr„_i 



dT\ I dk\...dki exp 



1 



~5>& 



l, in 



Here. 



I Im _ 'min[l,m] 



/ Ti Ti T\ 
T\ T 2 T 2 
T\ T 2 T 3 



M 

T 2 
T3 



V ri r 2 r 3 
To perform multiple integration over ki we use the well-known formula 



Tn J 



f +X l[dk 3 expl-i^fcjtO = 

J -°° 3 \ Z l,m J 



I (2tt)" 

drtT^"' 



It is straightforward to show that 

dctT< n > - Ti {t 2 - ri)(r 3 - r 2 )...(r n - r n _i), 
so that after integration over fej, Eq. I|27|l reduces to 

" ' A '' dr n _i n dr 2 p dn 



-9 

/27T 



<ir rl 



V 7 "" _ T n-1 -/0 V 7 ^ - 7"2 Jo V r l ( T 2 ~ T l ) 

=f ) / ^r:/ 2 - 1 !] / * m/9 - 1 (i-*r 1/9 «fe 

2tt/ Jo m=1 Jo 



(29) 



(30) 



(31) 



(32) 



(33) 



where g = gy/N/2. Thus, 



Z(g,N) = Z(g) = Yl 



(-.9)" 



_ e 9 



where 



'° r (I + ') 



*(«) = -F= / ' 



[!-$(£)], (34) 



(35) 



is the error function. 

When N — > oo, so does g\ expanding Z(g) for large g, 
we obtain 






Z(5) - ^ (---. + 0(r 5 )- 



Hence, 



71 = 1 + lim 



<91nZ 



at^oo 9 In N 



(36) 



(37) 



which is identical to the value we obtained in the previous 
section. 

Several remarks should be made at this point. First, 
note that while calculating the general term of the ex- 
pansion, we omitted the common factor 



Z (N)= / D[c(r)} e- H °M, 



(38) 



which actually has the form £ w (c.f. Eq. (jSJ). Sec- 
ondly, the coupling constant g plays here the role of the 
inverse cutoff length, which often occurs in field theory. 
Finally, we note that for any value of g > 0, we can 
find N large enough to make the non-dimensional cou- 
pling g = gwN/2 ^> 1, so that the number of acces- 
sible configurations scales as N" 1 ' 2 relative to the un- 
constrained case. This important observation is a signa- 
ture of universality: no matter how small g is, for long 
enough polymer the overall repulsion is infinitely strong! 



C. Probability distribution 

The Fourier-transform of the (unnormalizcd) proba- 
bility distribution function (PDF) for the end-to-end dis- 
tance of the random walk is given by 



G(d,N)= D[c(T)]e 



-iq-[c(J\T)-c(0)] e -H [c]-H I [fl,c] < 

(39) 

As above, we focus only on the z-dependent part of the 
PDF G(z, N). Expanding the path integral in powers of 
<?, and integrating out c(t) we observe that the ro-th term 
of the expansion reads 



n „n 

/ d"r.> 
2tt / J Jo 



dT n -\... I dn / Y\_ dk i ex P 



i 



n r N 



ar n / ar„_i... / -j== = exp 
o Jo 



o \/dctT( n ) 



I l,m 

~V-5>[TM] rm V ro ) 



<r_ 

2 



(40) 



where T^™^ is given by Eq. Ij32(l . Now, 



&l+i,- 



Si-i, 



[T{n)] -i 



n+i — n n — Ti-i 

' T 2 - Tl 
On—l,m 



Jm,l 



s u ( ? -4 Tr + yi 



0m,7i 



7~n—l T n T n — \ 



1 = 1 
I = n 



(41) 



Straightforward calculation yields 



Y,^ rin) V,r m =r n . 



(42) 



Equation l (4"0*|l therefore reduces to 



n r N 



dr„ 



dTn-l- 



dn 



VdetT( n ) 



exp 



-(JV- 



(-<?)' 



- g 2 AT/2 rN 



2"/ 2 r(n/2) 7 



r(n/2) 



MN/2)a 



(43) 



Summing over n, we obtain 



we finally obtain 



G( 9| JV) = e-^/ a ( 1 + 2 /* *«<«* W /($.)) ,(44) F (,, £) ^ A «) + -ML_ [°° *L (l - e - 2 ) 

\ Jo s / V27rg J Q u z \ J 



where 



Thus, 



/W = 4+xV (l-*(aO). 

\/7T 



-Z 2 /2N 

G(z,N) = —==-F&2), 



where z = z/yN and 

F(5,z) = l + 2 f 
Jo 



v^kN 



'^^^Vd-^l 



(45) 



(46) 



VT 



f(gs) (47) 



is the scaling function. To calculate F (g, z) we rewrite 
this expression as follows 

f 1 ds 1 - e-^'Vd- 3 ) 
F($,l)=A(s)-2 / == /($*), (48) 



where 



^(.9) 



1 + 2 



vr 



rf.s 



o sVT" 



/(fl«) 



(49) 



It seems that the integral in Eq. (|48|l cannot be calcu- 
lated analytically in general, i.e. for arbitrary g and z. 
However, since we are interested in the limit j> 1, we 
can easily calculate the leading term. We note that the 
integrand is essentially nonzero only for values of s larger 
than some value sq(z). However small Sq(z) is, we can 
always take g large enough to make gso(z) ^> 1. Thus, 
we can take 



f(gs) * - 



i 



2v / 7rgs' 



(50) 



so that 



, s Al s 1 f 1 ds 1 - e -l^ 2 /(i- 2 ) 
F(g,z)~A(g) + ^; == . (51) 



V^g Jo s 2 y/i _ s 2 



Using a substitution 



V2(l-« 2 )' 



(52) 



2irg J u z 

= A{g) 



V2g' 



For large values of g, we have 

1 



A{g) = 



2yRg 



;,2 



o(f 



(53) 



(54) 



Thus, when g — > oo, the scaling function is linear in z 
and the normalized PDF has the form 



G(z,N) = —e 

V ' ' 2N 



-z 2 /2N 



(55) 



Apart from a factor of 2, this function is identical to 
the one obtained in the previous sections. This factor 
appears because now the chain can be either in the z < 
or in the z > half-space. 




|z|/N' 



FIG. 2: The (normalized) scaling function F (g, z). Different 
curves are labeled by corresponding values of g. 

Figure shows the numerically computed scaling func- 
tion F (g, z) for different values of g. As expected, the 
larger g is, the closer is F (g, z) to the linear dependence. 



8 



IV. CONCLUSION 

To say the least, functional integration is not the 
most effective way to obtain the probability distribution 
G(z, N). Why then has one to work so hard if it is pos- 
sible to obtain the answer in just a few lines? 

Beside its clear educational value, path integral analy- 
sis of random walks is a much more flexible and powerful 
tool when it comes to real systems such as polymers in 
solvents. The self-avoidance constraint that we omit- 
ted so readily introduces long-range correlations that 
make the traditional approaches loose their power and el- 
egance. For example, factorization of G(r, N) into trans- 
verse and longitudinal parts is not valid anymore and 
integrating out one set of degrees of freedom introduces 
non-local interactions into the other. Whereas lattice 
walks are useful for Monte-Carlo simulations and dif- 
fusion equations could perhaps be modified to include 



mean-field corrections, field theory and primarily the 
renormalization group (RG) is today the only analytical 
tool to obtain universal scaling relations and phase di- 
agrams using controlled approximations 9 . For instance, 
generalizing our Hamiltonian to include self-interactions, 
interactions between a number of polymers, external po- 
tentials etc., it is possible to quantify the influence of 
space dimensionality and learn about the relevance of 
various interactionsii^. 
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